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We study the effects of a dark matter component that consists of bosonic particles with ultralight 
masses in the condensed state. We compare previous studies for both non-interacting condensates 
and with repulsive two-body terms and show consistency between the proposals. Furthermore, we 
explore the effects of rotation on a superfiuid dark matter condensate, assuming that a vortex lattice 
forms as seen in ultracold atomic gas experiments. The influence of such a lattice in virialization 
of gravitationally bound structures and on galactic rotation velocity curves is explored. With fine- 
tuning of the bosonic particle mass and the two-body repulsive interaction strength, we find that 
one can have sub-structure on rotation curves that resembles some observations in spiral galaxies. 
This occurs when the dark matter halo has an array of hollow cylinders. This can cause oscillatory 
behavior in the galactic rotation curves in similar fashion to the well-known effect of the spiral 
arms. We also consider how future experiments and numerical simulations with ultracold atomic 
gases could tell us more about such exotic dark matter proposals. 

PACS numbers: 95.35.+d, 98.62.Gq, 03.75.Lm 
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I. INTRODUCTION. 

In recent decades we have witnessed a breakthrough in 
precision cosmology with observational data from mea- 
surements of distant Supernovae [l|, Q and from WMAP 
3 indicating that our Universe is flat and dominated by 
non-standard forms of matter and energy. The data tell 
us that only about 4% of the energy resides in the bary- 
onic or luminous matter that we can observe through 
our telescopes. The rest of the energy is divided between 
the mysteries dark matter (26%) and dark energy (70%). 
The existence of dark matter (DM) was inferred earlier 
by its gravitational influence on the luminous matter in 
galaxies. Here the rotation velocity curves were found to 
disagree with expectations from general relativity. This 
lead people to introduce a dark matter halo around the 
galaxy that could explain the orbits. Models based on 
modification of the gravitational force are also being in- 
vestigated but thus far such an explanation seems im- 
probable Q. 

The nature of this non-baryonic DM is still a source 
of heavy theoretical and experimental investigation [5|. 
Over the years numerous matter particles have been sug- 
gested as candidates for the invisible dark matter. The 
most popular ones are the weakly interaction massive 
particles (WIMP) which have very low reaction cross sec- 
tions with normal baryonic matter Q , although, since it 
is non-zero, there are many proposals and on-going efforts 
to measure them directly [7[. One idea is that DM con- 
sists of scalar particles with extremely small masses, such 
as those appearing naturally as dilaton fields [81] in cos- 
mological applications of string/M-theory [9]. Ultralight 
DM was introduced in connection with a late-time cos- 
mological phase transition , and assumed to be a Bose 
liquid in the condensate state [ill ]. Later it was shown 
that problems with small-scale structure in the other- 
wise successful cold dark matter (CDM) model could be 
remedied by a Bose-Einstein condensate (BEC) DM com- 



ponent [12h15| |. since it will behave as non-relativistic 
CDM with a large quantum coherence length, avoiding 
the cusp behavior of DM halos found in simulations 13 . 
In Ref. [llj, and more recently in Refs. [16[ and 17j . 
galactic rotation curves were calculated for BEC DM and 
compared to observations of dwarf and spiral galaxies. 

Here we address the effects of rotation under the as- 
sumption that the dark matter halos are Bose-Einstein 
condensates (BEC) of ultralight particles. When a BEC 
is rotated at a rate exceeding some critical frequency, 
quantized vortices can be formed [l^ . The effects of 
vortex formation in galaxy halos has been addressed in 
Refs. [ill ■ In Ref. [ll| , the vortex density of the An- 
dromeda Galaxy M31 was considered, whereas Ref. [20| 
discusses a vortex as the cause of the flattened galaxy ro- 
tation curve. There has also been some interest in ana- 
logues between rotating spacetimes and superfiuid sys- 
tems [2l|. The latter is a continuation of recent claims 
by Chapline and collaborators that Black Hole event hori- 
zons represent quantum phase transitions as seen in Bose 
fluids and therefore should not be regarded as a break- 
down of classical general relativity |22l |. 



Vortex lattices are seen in experiments with BEC [23 
|26| when the sample is rotated at rates above the critical 
frequency. In those systems one has a superfiuid single- 
component BEC of alkali atoms at very low temperature 
with a two-body repulsive interaction described by the 
scattering length. Recently it was shown how an array 
of lasers can actually induce gravity-like interactions by 
avera ging over the usual electromagnetic dipole-dipole 
force [27| . This would in principle allow one to study 
the gravitational dynamics of BEC matter in the labo- 
ratory. Rotating such a system would allow us to study 
the properties of vortices when long-range 1 jr attraction 
is present and perhaps reveal how structures are formed. 
An alternative and likely more feasible way to approach 
the issue is through numerical simulations. Extensive 
theoretical and experimental work in ultracold atomic 
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gases has shown that Bose-Einstein condensates are accu- 
rately described by the non-linear Gross-Pitaevskii equa- 
tion [18| ■ Numerical simulations of systems governed by 
this equations including rotation and a gravitational po- 
tential would therefore presumably yield information on 
structure formation and vortex dynamics. 

In this article we want to explore some effects of rotat- 
ing dark matter halos that consist of ultralight bosonic 
particles in the condensed state. This sort of dark mat- 
ter would be superfluid and if rotated above the criti- 
cal angular velocity it would allow the formation of a 
vortex lattices in the halo cloud. Section |TT] considers 
the non-interacting case, whereas bosons with repulsive 
self-interactions are discussed in Section IIIII Effects of 
rotation on superfluid condensates are considered in Sec- 
tion IIVI in the light of previous works addressing this 
possibility in connection with halos of spiral galaxies. In 
Section [V] we explore the influence of vortex lattices in 
dark matter halos for early structure formation, in the 
virialization of gravitationally bound structures, and on 
the galactic rotation velocity curves. We compare the 
rotation curves in our ultralight BEC dark matter model 
to some observations of spiral galaxies and find that wig- 
gles in the rotation curves can be accomodated in our 
model within certain paramater regimes. An ultralight 
BEC dark matter component can thus be an additional 
source of oscillatory behavior in the rotation curves, al- 
though the main source is presumably the baryonic mat- 
ter in the spiral arms to which it is correlated [HI, H^] . 
Section IVT1 contains a discussion of some issues related to 
the fine-tuning of the mass of the bosonic particles and of 
the interaction strength of the repulsive two-body term. 
Some comparison between the presented astrophysical 
ideas and the current status of experiments and numeri- 
cal simulations with ultracold atomic gases are given in 
Section IVlTl Possible ways to test such exotic dark mat- 
ter proposals in the laboratory or on a computer are sug- 
gested. Finally, Section IVIIII contains a short summary 
and conclusions. 



II. ULTRALIGHT BEC DARK MATTER 

In this section we consider the possibility of a BEC of 
ultralight particles as the sole dark matter component. 
This was first suggested in Ref. [Tl| and further explored 
in Ref. [H| and Ref. Ql]. The idea is that ultralight 
particle to <C 1 eV will have very large de Broglie wave- 
lengths which means that quantum statistical effects are 
important and macroscopic coherent lumps of matter can 
emerge. As shown in Ref. [l9j], these light Bose particles 
will have a transition temperature to the condensed state 
that is of order T c ~ 2 • I0 9 K 1 , which is the expected 



This estimate of the critical temperature is based on comparing 
the thermal de Broglie wavelength [lSH to the mean inter particle 



temperature in the Universe after about I second. The 
fall of the temperature with the expansion means that an 
overwhelming majority of the Bose particles will be found 
in the condensate state. For the moment we will assume 
that the particles are non-interacting (we will relax this 
assumption below) and therefore only gravity acts on the 
system. Following Ref. [Hj], we now resort to Jeans insta- 
bility analysis to estimate our parameters. The growing 
mode under gravity is e 7 * with j 2 = AirGp, whereas the 
free field will be oscillatory; e~ lEt with E = k 2 /2m. The 
latter can be written as e 7 ' with j 2 = — (k 2 /2m) 2 . Not- 
ing that this is like normal Jeans analysis with sound 
speed c 2 = k/2m we have 7 2 = A-rrGp — (k 2 /2m) 2 . Set- 
ting this to zero, we get the Jeans scale 

rj = 2n/k J = TT 3 / i (Gp) 1 / 4 m-^ 2 

= 55m 22 1/2 (p/ (0fc )- 1 / 4 (fi m / l 2 )- 1 / 4 kpc, (I) 

where TO22 = to/10 -22 eV, and the background density 
is p b = 2.8 • 10 n ri m /i 2 M Q Mpc~ 3 . Below the Jeans scale 
the perturbations will be stable and above it they will 
behave as ordinary CDM [l^]. The stability below the 
Jeans scale is guaranteed by the uncertainty principle. 
If the particles are confined further, their momenta will 
increase and oppose the gravitational contraction. 

From the Jeans scale in Eq. (fTJ) we can see that the 
mass has to be extremely small for the current scenario 
to be responsible for galaxy-size structures at present. A 
typical galaxy has a matter density of about 200p(, and a 
radius of order r ~ 10 kpc. Using = 0.3 and h = 0.7 in 
Eq. |T]), we see that a mass of precisely order m ~ 10~ 22 
eV is needed for the condensed particles. As noted in 
Ref. [Hj], the observational evidence can accommodate 
even lighter fields m < 10~ 33 eV in quintessence models. 
Ref. [12| proceeds with one-dimensional simulations and 
show that such ultralight Bose particles could actually 
help solve cusp problem in dark matter halos by sup- 
pressing small-scale linear. This was confirmed in more 
detail in Ref. [30|. 

As mentioned in the introduction, the rotation curve 
of galaxies is an important feature that dark matter ex- 
plains very well. In Ref. [l9j] the ultralight BEC dark 
matter scenario was used to predict the rotation curve 
and the results were compared to observations from 
the Andromeda galaxy M31. The model was found to 
agree with the data if the boson mass was in the range 
to ~ 10~ 24 — 10 -23 eV, in rough agreement with the 
analysis in Ref. [l2l ]. 



distance. For condensation to set in these two quantities must 
be similar, giving the relation ksT c ~ fecfp/m) 1 / 3 where p is 
the mean density of the Bosonic matter [19(1 . A more detailed 
calculation of the transition temperature for a uniform Bose gas 
yields the expression ksT c ~ /i 2 (p/m) 2 / 3 /m. Using the latter 
relation yields an even higher value for T c and does not alter the 
argument in the text that the Bose particles are expected to be 
in the condensed state at the present epoch. 
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III. INCLUDING SELF-INTERACTION 

There has been a number of investigations into dark 
matter models with Bose particles that consider the 
(likely) more realistic case of bosons with repulsive self- 
interactions [HI [lij]. The motivation for this was the 
problem of too much sub-galactic structure that was 
mentioned above. The typical case considered is that 
of a quartic self-interaction A</> 4 (A > 0) which is the 
generic textbook example of an interacting field theory. 
From an atomic physics point of view the quartic case 
is also very interesting as the mean-field approximation 
yields the Gross-Pitaevskii equation (GPE) for the con- 
densate state that has proven very useful in describing 
BEC experiments [l8j . 

Even more interesting is the fact that there is actually 
an exact solution for the case of a self-gravitating Bose 
gas with a repulsive quartic term (l3l Uj\ . This comes 
about since in this case the GPE can be recast into the 
Lane-Emden equation with polytropic index n = 1 which 
has the analytic solution 

.W^o^^v^G. (2) 

r/r 

K comes from the equation of state: p = Kp 2 . This 
latter relation can be found by simply estimating the 
ground-state energy of a Bose condensate with repul- 
sive self-interaction 18]. In atomic physics, one would 
typically relate the interaction strength (specified by A 
above) to the scattering length for two-body scatter- 
ing at very low energy. In our current notation this 
yields K = gn/2m 2 , n being the particle density and 
g = Aith 2 a/m is called the interaction coupling constant 
pL8j ] . The sound velocity in the gas is given by c 2 = gn/m. 
The observed dark matter halos do not have the profile 
given by Eq. J2j but, as noted in Ref. [ljj, there could be 
non-zero momentum particles present that would allow a 
more realistic power-law behavior outside the core. 

In this repulsive self-interaction scenario one can es- 
timate the mass of the Bose particles by using the first 
zero of Eq. ([5} , R — nr , which is assumed to be the halo 
radius. Combining all the formula given above yields 



GR 2 



1/3 



= 6.73- KrV /3 i?- 2/3 eV, 



(3) 



where a is measured in fermi and R in kpc. Taking R ~ 
10 kpc and a between 1 fm and 1 nm (typical atomic 
magnitudes fl8T]). we obtain masses in the range from 
m ~ 10 meV to 1 eV. We immediately notice that this is 
many orders of magnitude away from the non-interacting 
case of the last section. A mass of m ~ 10 -22 eV would 
require an extremely low scattering length of a ~ 3.3 • 
10~ 61 fm (keeping the same halo size). This is of course 
extremely small and borders on the non-interacting case 
which tells us that there is consistency with the results 
presented in the previous section. 



At this point one might wonder if the critical temper- 
ature estimates given earlier will be altered significantly 
by the fact that the bosons now have a repulsive two- 
body interaction. In the context of the uniform Bose gas 
model, this has been addressed in Ref. [3l| where is was 
shown that the change in the critical temperature obeys 



6T C /T C 



an 



1/3 



If we measure a in fermi, m in eV, and 



p in 10 23 g/cm 3 (typical halo density) then we get 
8T C /T C ~ 2 • 10- 10 a(p/m) 1 / 3 . We therefore see that the 
correction to T c is extremely small for both m ~ 10~ 2 — 1 
eV and for ultralight m ~ 10 -22 eV (remembering that 
a oc m 3 according to Eq. ([3])). We therefore see that the 
estimate T c > 2 • 10 9 K given earlier will not be altered 
by the repulsive self-interaction. 

The considerations in this and the previous section 
rely on Newtonian gravity, and general relativistic ef- 
fects could modify the results. A general relativistic de- 
scription of the BEC dark matter scenario is discussed 
in Ref. 17J, particularly the effects on the rotational 
curves of galaxies. It was found that rotational curves 
of a number of galaxies could be well reproduced, al- 
though the Newtonian analysis already produces a tan- 
gential velocity of test particles at the halo radius of 365 
km/s, consistent with observations. Note also that we 
are tacitly assuming that the ultralight Bose particles 
are non-relativistic which was shown to be a very good 
approximation in Ref. [ill . Il9j . We therefore restrict our 
discussion to Newtonian gravity here. 



IV. THE EFFECTS OF ROTATION 

Several studies have addressed the effects of rotation 
on BEC dark matter. In Ref. [ijj there is a brief discus- 
sion of the effect on the Lane-Emden equation, whereas 
Ref. [2(| consider a BEC of axions with a single vortex 
arising from global rotation in the early Universe. This 
latter scenario is, however, less likely to occur since the 
global rotation rate of the Universe can be estimated from 
various observations and is likely very small if non-zero 
[38| . The formation of vortices in condensates is known 
from atomic physics [l|[. For superfluid Bose systems 
any angular momentum imparted will reside in quan- 
tized vortices that constitute small regions where the su- 
perfluid density goes to zero on a scale comparable to 
the coherence length (to be discussed below). As already 
noted in Ref. [l3j BEC dark matter with self-interactions 
will actually constitute a superfluid. 

A superfluid will only form these vortices when the 
rotation rate exceed a critical frequency £! c . For a sample 
of size R with coherence length £ the critical value is [l8l ] 



mR 2 



In 



1.46i?\ 



6.21 -10" 17 , flA6R 
In 



m 2 2 



R 2 



( 4 ) 

where 77122 = m/10 22 eV and R is measured in kpc. 
This critical frequency is based on the Gross-Pitaevskii 
equation and the coherence length £ is given by the self- 
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interaction through £ = l/y/8iran with n the particle 
density and a the scattering length. The formation of 
vortices has been experimentally demonstrated in super- 
fluid 4 He [23[. Later they were produced in dilute BEC 
alkali atom gases f2~4T - |26j where beautiful lattices where 
created with more than 100 vortices in regular align- 
ment. These experiments also showed a density profile 
that had dips at the vortex positions, confirming the the- 
oretical expectation of density drops in the vortex core 
mentioned above, although the actual experiments had 
cores with small but non-zero density. The expected in- 
crease in the number of vortices with increasing rotation 
frequency was also observed. Theoretically, the vortex 
line density is given by n v — mfl/Tth which yields a total 
number of vortices N v — irR 2 n v . 

The idea of similar vortex formation in rotated BEC 
dark matter was first considered in Ref. Here the 

rotation rate of the Andromeda galaxy M31 was calcu- 
lated and compared to the critical frequency. M31 was 
found to rotate at a rate vastly above fi c and vortices 
were therefore suggested as a possible consequence. The 
total number of vortices for M31 was found to be around 
N v ~ 340. As a possible means of detecting these vor- 
tices, Ref. [l9[ suggested gravitational lensing or polar- 
ization effects from frame-dragging of light from distant 
background sources. The use of gravitational lensing to 
infer the possible existence of a BEC dark matter com- 
ponent was also suggested and explored in Ref. [l7j . 

As already mentioned above, the halo radius can 
be taken as the first zero of the solution in Eq. ©, 
i.e. R 2 = TT 2 h 2 a/Gm 3 . In the scalar field treatment 
with interaction A</> 4 , the corresponding size estimate is 
R 2 = TT 2 \h 3 /8irGm 4 c [lji |32[. Equating these, we ob- 
tain A = Snamc/h = 1.3- 10" 29 (a/l fm)(m/10- 22 eV). 
Using m — 10~ 22 eV, a halo mass density phaio ~ 10~ 23 
g/cm -3 , and assuming that £ ~ 1 kpc, we find A ~ 10 -92 , 
a very small value. As noted earlier, the repulsive self- 
interaction is not likely to influence the gravitational for- 
mation of the halo 15]. However, as discussed above 
such small values give large vortex cores when fi 3> f2 c . 
In Ref. [32j j . fits to rotation curves were presented and 
it was found that m 4 /X - 50 - 75 eV 4 . For m = 10~ 22 
eV this gives A ~ 10 -90 , close to the present estimate. 
Importantly, our smaller A actually fulfills the Big Bang 
Nucleosynthesis (BBN) bounds discussed in Ref. |32j ]. 



V. A DARK MATTER VORTEX LATTICE 

If dark matter contains a component of condensed 
BEC particles that is superfluid and if the halos are ro- 
tating then it is not inconceivable that there can be vor- 
tex formation as discussed above. However, the quan- 
tized vortex discussion of Ref. [19( makes an important 
assumption about the coherence length, £, entering Q c 
in Eq. (jj). £ is taken to be of kpc size. This is based 
on the non-interaction arguments of Section [TTJ With 
no self-interaction there is only the gravitational scale 



h 2 /GMm 2 available, which becomes of galactic size for 
masses m ~ 10~ 22 eV. However, when including self- 
interactions as in Section Hill through a scattering length 
a, there is also a scale given by £ = l/%/87rar7, which is 
the usual Gross-Pitaevskii coherence length. The latter 
coincides with the characteristic length over which the 
density is expect to go to zero in a vortex. 

In the following discussion we will make the assump- 
tion that it is actually the coherence length £ that deter- 
mines the vortex core size. This is reasonable since the 
vortices are local entities and as such should depend on 
the local interactions. From this point of view the gravi- 
tational scale H 2 /GMm 2 merely serves to determine the 
total halo cloud size. The coherence length can be writ- 
ten as 



£ = 2.73 • IQT^mga-WpV? kpc, (5) 



where pdm is the dark matter halo mass density mea- 
sured in 10~ 23 g/cm 3 , a is in Fermi, and 77122 = m/10~ 22 
eV. The measure for pdm originates from assuming that 
a halo has a radius of about 10 kpc and contains about 
3-10 11 M of mass giving a density of about 1.74-10 
g/cm 3 . The front factor in Eq. (0) seems exceedingly 
small and it would seem that any vortices on this scale 
would be completely irrelevant to the galaxy structure. 
However, we recall from the mass estimate in Eq. ([3]) that 
a also needs to be extremely small. Inserting a value of 
a ~ 10~ 61 fm (corresponding tom~ 10~ 22 eV) actually 
gives £ ~ 0.05 kpc, so that the vortex size is sub galac- 
tic but sizable. This demonstrates that vortex structures 
could be important for the structure of rotating galaxies 
and we will discuss possible effects below. 

Before we consider vortex lattice effects on structure 
and evolution, we now briefly compare our suggestions 
with the observational data from M31. Ref. [l9| estimate 
the critical frequency to be about f2 c = 2 ■ 10~ 19 rad/s 
for M31 (using R ~ 150 kpc and m ~ 2 • 10~ 23 eV). 
The estimate for the actual rotation frequency of M31 
based on observations is given as f2 = 5 • 10~ 17 rad/s, 
so that the assumption O 3> f2 c is justified and a vortex 
lattice is therefore possible. However, in Ref. [lj| the co- 
herence length used in Eq. Q was taken to be £ ~ 30 
kpc (the gravitational length scale for the non-interacting 
condensate discussed above). We are assuming that the 
coherence length scale for vortex formation is that as- 
sociated with the two-body repulsion. This is, however, 
not a serious obstruction since it enters only in the log- 
arithmic term in Eq. Using the value quoted above 
of £ ~ 0.05 kpc in the critical frequency would there- 
fore only amount to about a factor of five increase which 
would not jeopardize the relation Q 3> f2 c . Furthermore, 
changing the two-body repulsion (through changes in a) 
would not be severe for the same reason. 
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A. Influence on Early Structure Formation 

For the vortex lattice to appear in the current scenario 
with BEC dark matter we need rotation. One can as- 
sume that rotation is a primordial feature of the Universe. 
This possibility has been explored by many authors (33l — 
. However, this is at odds with considerations of CMB 
anisotropies [37|, which sug gest that the primordial rota- 
tion rate is very small [38|. The rotation that one finds 
in galaxies is therefore not considered to be primordial 
and is believed to arise from gravitational tidal torque 
forces during the growth of non-rotating initial pertur- 
bations toward virialized galactic structures. This way 
of generating angular momentum is well studied [39-41] 
and widely accepted as a key mechanism. 

Although primordial rotation seems out of the ques- 
tion, we will briefly consider what would be the result if 
it was found that our Universe could have been rotating 
in early epochs. More precisely we want to explore what 
would be the effect of a rotation rate large enough to 
support a vortex lattice at the time of decoupling when 
baryonic perturbations can start to grow. We can ar- 
gue in the usual way that when baryonic perturbations 
can finally grow, they will do so under the influence of a 
potential well that is created by the dark matter pertur- 
bations. Assuming that the baryons themselves provide 
a negligible contribution to the background density, a 
Jeans analysis tells us that the baryonic density contrast 
in the linear regime obeys [42j 



we assume that these can (at least locally) be neglected. 
This assumption would of course have to be checked in 
more detailed investigations. 

Since there is no evidence for primordial rotation, the 
discussion above is likely of little relevance to the Uni- 
verse we live in. The rotation that we observe in spirals 
must therefore originate in late-time events when the per- 
turbations have entered the non-linear regime and Jeans 
analysis is no longer useful. Estimates show that the 
galactic structures that we see today cannot not have 
formed at redshifts much larger than z ~ 10 43]. Lin- 
ear perturbation theory is expected to hold around de- 
coupling which is z ~ 1000, so we are well beyond this 
approximation. The structure of baryons is also com- 
plicated by pressure and radiative terms, which means 
that hydrodynamical codes must be used. State-of-the- 
art in structure formation are highly involved iV-body 
computer simulations that can calculate the evolution of 
the initial spectrum of perturbations into the non-linear 
regime. However, these simulations do show that the 
gravitationally bound systems that arise tend to virial- 
ize on a fairly short timescale. In this article we will not 
consider such advanced simulation techniques but merely 
consider what the effects of vortices could be on the viri- 
alized structures, since these are the ones that we observe 
to be rotating. 



Influence on Virialization 



■i 4 i 



k B T \ k^ &B = 25dm (g) 
m p J a 3 3 t 2 ' 



where t is cosmic time, a is the scale factor, T the temper- 
ature of the baryons, m p the proton mass, kg the Boltz- 
mann constant, and k the wave vector of the Fourier den- 
sity mode under consideration. For a matter-dominated 
Universe with 5dm °^ a i one finds 



S B (t) 



'DM 



(*) 



1 + Ak 2 ' 



(7) 



where A = 3/2(kBT /m p )(t 2 /a 3 ) with T the tempera- 
ture today (using matter-dominance from decoupling to 
the present). Baryons are thus coupled tightly to dark- 
matter on large scales, whereas the pressure support on 
small scales (Ak 2 <C 1) suppresses growth. If there was a 
primordial rotation of the dark matter of sufficient mag- 
nitude to cause vortex formation, then one would have to 
modify the driving term 5dm on the right-hand side of 
Eq. to reflect the vortex structure. Given some char- 
acteristic vortex size (as in Eq. ([5])) one would have to 
make the Fourier transform of a swiss-cheese like config- 
uration of the dark-matter. This would then give a new 
driving term and thus influence the structure formation 
in the baryons. There is a technical point here that we 
have ignored, which is that the Jeans analysis does not 
strictly apply to rotating systems. There would be cor- 
rections to this from Coriolis and centrifugal terms but 



Let us now address the question of the virialization 
process itself. The usual argument is that the bound 
system collapses under gravity with internal heating as 
a result. The virial theorem tells us that the internal 
kinetic energy should be half the gravitational potential. 
The systems therefore collapses until a radius is reached 
where this condition is fulfilled. However, when we in- 
troduce short-range two-body repulsion for the BEC dark 
matter, the virial theorem is modified and becomes 



2E, 



kin 



E, 



grav 



3E re p 0, 



(8) 



where E^m is the kinetic energy, E grav the gravitational 
potential energy, and E rep is the energy arising from the 
two-body term. The latter contribution is easy to calcu- 
late for our zero-range potential and becomes 



E re p — 



2nh 2 a 



dfn 



(9) 



where n(f) is the density of BEC dark matter. ./V-body 
simulations and observations tell us that objects virial- 
ize. So the contributions from E rep must be negligible 
in virialization to not jeopardize this fact. Consider the 
ratio 



E, 



rep 



E„ 



(10) 



Here we are assuming a uniform sphere for the matter 
distribution which gives E grav = —3/5GM 2 /R. In order 
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for E rep to not influence virialization, we need this ratio 
to be much smaller than 1. However, since we have the 
relation between a and m in Eq. (j3|) this ratio becomes 
simply — 5/27T 2 ~ —0.25, so we see that the repulsive 
interaction will be non- negligible. We therefore have to 
reconsider how to apply the virial theorem. This we will 
do along the lines of the simple top-hat model as dis- 
cussed in Ref. [43| . 

The virial theorem with inclusion of the repulsive two- 
body interaction and the use of Eq. for m yields 



Ekin — {~\~Eg rav 3E re p) — 



1 

2 

-0. \ c 2iEg Tav . 



2 An 2 J * grav 
(11) 



We see that the internal kinetic energy is a consider- 
ably smaller fraction of the gravitational energy than for 
standard virial consideration with a = 0. If we have a 
sphere of matter with radius r max that starts to collapse 
from rest, then from conservation of energy one can eas- 
ily estimate that the virial theorem will be satisfied when 
T ~ r maa; /1.14, whereas for a = is would be r = r max /2 
^43] . Virialization thus predicts a less compact system in 
the a 7^ case. This means that when gravitational ob- 
jects become virialized they will have densities that are 
about 27pb, as compared to the a — estimate of 150pb 
[43j . The two-body repulsive interactions are thus seen 
to yield less dense objects, something that could be very 
good in terms of getting less dense halos than some sim- 
ulations produce and that are at odds with observations. 
To test whether this helps would require full simulations 
with ultralight BEC dark matter. 

iV-body simulations predict that virialization occurs 
at densities above 150pt, (estimates put the value close to 
400/Ot, [H|])> which is an order of magnitude larger than 
the calculation above. Ultralight BEC dark matter would 
therefore seem to be ruled out. However, the above re- 
sults are only as good as the assumptions we use in their 
derivation. In particular, the use of Eq. is based on 
the Lane-Emden result of Eq. ©. We know that halos 
do not have this profile. We should therefore not put to 
much emphasis on the mass m in Eq. ([3]) , but rather leave 



m as a parameter. If we assume that E r 



-kE„ 



one can invert the virialization analysis of Ref. [43j of to 
get an estimate of k if we are to get virialized densities 
of > 400p 6 . This yields k < 0.1. For this to be fulfilled, 
Eq. dTU]) tells us that ajm\ 2 < 4 • 10~ 44 fm/eV 3 with 
a in fermi and 77122 = m/10~ 22 eV. This shows that we 
can actually accommodate TV-body results on virializa- 
tion are fulfilled by having a small value of a for ultralight 
Bose particles. Such small values for a are also consistent 
with the values predicted from the arguments in earlier 
sections of a ~ 10 -61 fm. 



C. Influence on Galactic Rotation Curves 

Another interesting issue is whether a dark matter 
vortex lattice can influence the galactic rotation curve. 



These curves are key pieces of evidence for the existence 
of dark matter due to their flat (or slightly increasing) 
behavior at large radii where standard Newtonian theory 
including only luminous matter predicts a sharp decrease 
with radius. The galactic rotation velocity is calculated 
from the basic formula 42] 



v 2 (r) 



GM(r) AttG 



dxx 2 p{x), (12) 



where p(r) is the total mass contribution from dark and 
luminous matter. In the present work we assume that 
the dark halo contribution has the form 

Po 



Ph{r) 







1 + 


fey 







(13) 



where r$ is related to the core radius and 7 is the expo- 
nent determining the large distance behavior. This form 
avoids the small distance cusp that BEC DM presumably 
solves (l2T - fl5j . The rotation curve can now be calculated 
from the simple formula v 2 — v 2 m . k + vf um , where i>i um is 
the contribution from the luminous matter in the bulge, 
disk, and surrounding gas. 

Experiments with rotating atomic BEC produce tri- 
angular vortex lattices. Theoretically, this is is the 
minimum energy configuration, with the square lattice 
slightly higher in energy (44j. Recent studies including 
long-range dipole-dipole interactions between the bosons 
find that other configurations are preferred as the long- 
range forces increase po| . and we speculate that similar 
results can arise for the 1/r gravitational force. How- 
ever, since it is unclear which kind of lattice would arise 
in the galactic halos, we have explored both triangular 
and square lattices. We also mention that square lattices 
are commonly seen in two-component condensates [46| . 
which would be relevant for multiple species BEC DM. 

The vortex lattice is modeled as an array of zero- 
density tubular regions of square cross section in the DM 
mass density. We have thus ignored the spherical shape 
and the smooth fall-off to zero density that the more re- 
alistic vortex cores would have. We have checked that 
this simplification did not have significant effects on our 
results. Due to this implementation, the vortex size rep- 
resents an effective coherence length, simply related to 
£. Keep in mind that the square shape can cause some 
rather sharp features in the calculated rotation curves 
which are smoothed for more realistic profiles. To find 
the rotation velocity we have numerically integrated this 
'swiss-cheese' halo along with the luminous contributions. 
In effect the method used assumes that the lattice rotates 
with the system, which is what is seen in cold gas exper- 
iments as well. To demonstrate the effect of the vortex 
lattice, we show in Figure []] a number of rotation curves 
based solely on the DM. These have been calculated for 
vortex size 5 kpc and vortex-vortex distances of 1 and 
2.5 kpc. These parameters are optimal in order to get 
the oscillations at short distances shown in Figure [TJ For 
smaller core sizes the effect vanishes and for larger ones 
the scale of wiggles is too large to fit observations. 
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FIG. 1: Rotation curves based on BEC DM with a vortex lattice for different halo sizes ro = 8 (full), ro = 4 (dotted), and 
ro — 2 kpc (dashed). The vortex core size is 5 kpc. In the left column the distance between vortices is 2.5 kpc, whereas it is f 
kpc in the right column. Upper row is for a triangular lattice and the lower one is for a square lattice. 



At this point we mention that we have used a cutoff on 
the vortex lattice for small radii (of order the halo size 
ro), effectively ignoring vortices that are located in the 
halo center. With this cutoff we can fit the observation 
of large constant density in the center of the galaxies. 
This appears reasonable since luminous matter should be 
dominant in the inner-most region. Another issue is the 
density of vortices. In our simple model we have assumed 
this to be constant with the use of regular lattices but 
there could be variations [47j. 

The wiggles found in the velocity curves with vortices 
present are quite interesting in relation to observations. 
Here one also sees distinct wiggling features in the plots 
[U G3- In Ref. [H the data from NGC2998 was ex- 
plained by postulating the BEC dark matter to be in an 
excited condensate state which has several zeros in the 
wave function. When this is translated into a dark mat- 
ter distribution it gives areas of low density in the halo. 
In the velocity curves this will of course have similar ef- 
fects as the vortex lattice with its empty cores. However, 
with vortices one can get wiggles without having to ex- 



plain how the excitation occurred and why the conden- 
sate should be in a particular excited state at present. 
One would of course also have to include the luminous 
matter for a real comparison to data, and for small radii 
this will presumably dominate the rotation curves [48j . 
Below we will include the effects of the luminous compo- 
nent explicitly when comparing to observations. 

Before embarking on a comparison to observations, we 
note a number of problems with this model. First, we 
have assumed that the vortex cores have the same shape 
at all radii, independent of the local halo distribution. 
This is likely unrealistic as we would expect some effects 
of the finite size toward the edge. Secondly, for the co- 
herence length to approach 5-10 kpc, from Eq. f5|), we 
would need either a ~ 10 -65 fm or an increase in 77122 by 
four orders of magnitude (or a combination of both). We 
estimated earlier that existing suggestions for the param- 
eters of ultralight BEC dark matter of a ~ 1CP 61 fm and 
m ~ 10~ 22 eV gives £ ~ 0.05 kpc, which would be im- 
possible to see in velocity curves. So it would seem that 
we need some fine-tuning of m and a to get any effect. 
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This we will address below. 



D. Comparison to Spiral Galaxies 

An interesting feature of the BEC DM model is the 
oscillations in the rotation velocity at small radii, and we 
now compare the model to spirals that show pronounced 
oscillatory features. Both dark and luminous matter play 
important roles in modeling spirals, and we therefore in- 
clude the luminous component based on recommended 
values from various observation as described below. Note 
that we allow the dark matter component and the details 
of the vortices to vary in order to best describe observa- 
tions. For simplicity, we ignore the feedback effect on the 
luminous matter of changing the dark component. Since 
luminous matter dominate the rotation curve at small 
radii and dark matter at larger radii in our model, this 
should be a fair approximation. 

We consider a sample of spirals that have been used in 
BEC dark matter studies previously [lj [H, frit. I n Fig- 
ure [2] we show data for the Sc spiral NGC2998 The 
bulge was assumed to give a constant contribution for 
r < 3 kpc and then fall as r -1 / 2 whereas the disk and 
gas are combined into the functional form (r/ (r + b)) 6 of 
range b [13] • The magnitudes of these non-dark contribu- 
tions was scaled to reproduce those in Ref. [H(J, whereas 
the DM contribution is varied to best describe the ob- 
servations. It can be seen in Figure [5] that a square lat- 
tice is closer to the data. In particular, we see good 
agreement between maxima and minima of the curves, 
although the magnitudes are not well reproduced, par- 
ticularly at larger radii. We find that the oscillations in 
the model are directly related to the vortex core size, L v , 
and lattice distance, d, and the characteristic length of 
the wiggles observed in the rotation curves can thus be 
matched in the present model by careful choice of pa- 
rameters. However, as explained before the main source 
of the oscillatory behavior is considered to be baryonic 
through due to the presence of spiral arms. 

To further explore the model, we consider observations 
from two Sc spirals, NGC753 and NGC801 (H, a dwarf 
smral, NGC1560 (HJ], and a thin disk spiral, NGC3198 
[52j in Figure [3] The non-dark components are mod- 
eled as for NGC2998 above. Notice that NGC1560 and 
NGC3198 have no bulge. The DM for NGC753 and 
NGC801 assumes a square vortex lattice (a triangular 
lattice compares worse with data), and we see that this 
gives a better agreement with the observed wiggles than 
for NGC2998, but still with discrepancies between min- 
ima and maxima. This could likely be resolved, if one 
makes the vortex density non-uniform as suggested in 
Ref. [47j. In the dwarf case of NGC1560, we see that 
our model can reproduce the observed kink in the curve, 
in this case even better with the triangular lattice. This 
feature was not reproduced in previous fits to BEC DM 
[HE3- In the thin disk s P iral NGC3198 our model only 
gives a small improvement over previous studies, and 



square and triangular results are practically the same. 
There are some discrepancies at large radii, which would 
likely be reconciled by a better description of the gas 
component. 

The BEC DM model can accomodate data on differ- 
ent spiral types, and the parameters used are very simi- 
lar. The vortex cores are 5.0-6.5 kpc in size with a lattice 
spacing of 1-2 kpc. We note that experiments have shown 
that the core and vortex- vortex distance can be of similar 
magnitude [53| in trapped condensates, although the core 
size is smaller than the distance. The theoretical vor- 
tex density expected in uniform condensates is given by 
n v = mfl/irti = 0.05(r7i/10- 22 eV)(f7/10- 16 s- 1 ) kpc~ 2 . 
For the spirals considered here we get n v ~ 0.10 — 0.13 
kpc -2 . The total number of vortices is then N — 
wR 2 n v ~ 30 - 1000 for radii R ~ 10 - 50 kpc. For 
the square lattice this gives about 45 vortices within a 
distance of 50 kpc using the size and distance quoted 
above, which is a reasonable halo size for the Sc spirals 
discussed. Our model parameters are thus not completely 
untenable. With smaller cores we could get larger den- 
sities or have larger distance between vortices. However, 
we have found that in order to obtain the wiggles in the 
rotation curves the core size has to be roughly equal to 
the length between maxima and minima in the observa- 
tions. In NGC3198 we have used a larger distance of 2.0 
kpc to get a good fit, consistent with its rotational ve- 
locity being slightly lower than that of the others. Also 
for the dwarf NGC1560 we have used a smaller core size. 
This gives a larger vortex density which is consistent with 
a slightly larger rotational rate compared to the Sc sam- 
ples. 

The wiggles in the rotation curves are located at short 
radii and studies have shown them to be correlated with 
the spiral arms [H, [Hj], both facts strongly indicating 
baryonic processes as the origin. The current model im- 
plies that there could be additional effects in an ultralight 
BEC dark matter scenario. Assuming that the wiggles 
are caused by baryonic matter only, we can use the cur- 
rent model to extract bounds on the parameters from the 
lack of vortex lattice effects in observations. This would 
in turn influence other models of BEC DM through lim- 
its on boson masses and two-body coupling terms. In 
particular, the absence of large wiggles produced by vor- 
tex lattices would give a lower limit on the two-body 
repulsion of the bosons. 



VI. FINE-TUNING ISSUES 

The parameters of importance to the success of an 
ultralight BEC dark matter model are naturally the 
mass, which must be extremely small, and, for the self- 
interacting scenario, the strength of the two-body repul- 
sion between the condensed particles. In Ref. [l9] BEC 
dark matter is suggested to arise from a single scaler 
field coupled to gravity undergoing spontaneous symme- 
try breaking to acquire a vacuum expectation value. This 
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FIG. 2: Rotation curve data for the Spiral galaxy NGC2998 (dots with uncertainties) and calculation with a square (thick line) 
and triangular (short dashed) BEC vortex lattice in the DM component. Shown are also the individual contributions from the 
bulge (dot-dashed), disk+gas (dotted), and DM for the square lattice case (long dashed). The parameters used are ro = 5.0, 
7 = 2.15, b — 3.0, L v — 6.5, and d — 1.0 (all lengths in kpc). 



produces a cosmological constant A. The breaking of 
symmetry is done via a Ginzburg-Landau potential with 
quadratic and quartic terms. This gives mass and inter- 
action terms to the scalar field. If we now make the ad- 
ditional assumption that the vacuum expectation value, 
</>o, arises from a mechanism that preserves parity (so 
that third-order terms can be ignored), the interaction 
term is 



(mc 2 ) 2 ^ 4 _ <j/_ 
<±{hcfcg ~ 5 T' 



(14) 



where we have restored constants of c and h. This 
terms is of course merely the standard interaction term 
in the Gross-Pitaevskii theory of interacting condensed 
bosons. We therefore see that the self- interacting sce- 
nario emerges from this procedure. 

The above scenario actually contains some additional 
information about the mass and interaction strength of 
the scalar field. One can derive a relation between mass, 



the interaction strength, and the cosmological constant. 
As in the previous sections, we want to use g = 4jvh 2 a/m 
to express the interaction through the scattering length 
a. We have 



4tt^ 6 A 



= 7-10 



A 



(15) 



where A = 0.7p cnt 
the cosmological constant from WMAP Q . As is clearly 
seen, the value of a in this model is quite far from the 
considerations of the previous sections, where a ~ 10~ 61 
fm was found to be a favorable value. A mass of order 
m ~ 10 -28 eV would be required to get this latter value 
of a (this is still within the bounds on ultralight dark 
matter discussed in Ref. (EH). 

These considerations imply that the dark matter vor- 
tex lattice proposed here is far from generic and requires 
finely tuned parameters to work. However, the model 
of spontaneous symmetry-breaking discussed in Ref. [l9T | 
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FIG. 3: Rotation curves for NGC 1560 (5.0,2.15,0.5,5.0,1-0), NGC753 (5.0,2.00,3.0,6.5,1.0), NGC801 (5.0,2.15,3.0,6.5,1.0), and 
NGC3198 (5.0,2.00,1.2,6.5,2.0) compared to model calculations. The plot assignments are as in Figure(2] The parenthesis give 
the parameters used in the format (ro, 7, b, L v , d)). Notice the change of radial distance scale for NGC1560. 



is also highly speculative. The cause of the observed 
dark energy is still unknown and it is therefore reason- 
able to consider scenarios where the relation in Eq. (fT5|) 
is not valid. This can easily be done by considering more 
general potentials than the standard Ginzburg-Landau 
scalar potential. The breaking of symmetry could then 
be accomplished in a manner that allows other relations 
between m and a. This could allow other ranges of the 
parameters and perhaps make way for dark matter BEC 
vortex lattices without extreme fine-tuning. For lack of 
a better understanding, the cosmological constant, A, 
could be considered an outside parameter, we want to 
address dark matter scenarios that are not necessarily 
from unified dark matter-dark energy models. 



VII. RELEVANT EXPERIMENTS AND 
NUMERICAL SIMULATIONS 

The present experimental situation in ultracold atomic 
gases also encourages the hope that one could possibly 
probe the consequences of bosonic dark matter and con- 
densates in laboratory experiments. Many groups around 
the world can routinely produce condensates of bosons, 
degenerate fermi gases, and also interesting mixtures of 
the two [H, Hfil • Another feature is the extreme control 
experimentalists exercise over the interactions between 
the different atomic species. This is achieved through 
the use of Feshbach resonances that allows the tuning of 
the interactions strength over many orders of magnitude 
and also whether the atoms repel or attract each other. 

A very interesting proposal is that one can tune the in- 
ter atomic dipole-dipole interaction to resemble the grav- 
itational force [27]. By illuminating the atomic clouds 
with a careful arrangement of multiple laser beams, it 
was suggested how to eliminate the 1/r 3 dependence of 
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the dipole-dipole interaction and keep only the 1/r part 
in the near-zone. The force constant of the remaining 
term depends on the laser intensity and atomic polariz- 
ability. The resulting force is equivalent to the attraction 
between two opposite electric charges with q ~ e/2000, 
so compared to electromagnetism it is not strong. How- 
ever, in comparison, the tiny magnitude of gravity means 
that this force can be much larger than the normal grav- 
itational force based on the atomic masses. 

If this proposal is successfully implemented in exper- 
iment, one could therefore simulate gravitational forces 
between the gas particles of both fermi and boson species, 
and with 'effective' masses that are vastly different from 
the given atomic masses. This could very likely be a 
way to test the movement of particles during virializa- 
tion and see how structures form. Comparison with ob- 
servations and iV-body simulations would then allow us 
to expand our knowledge of the evolution of structure 
in our Universe. Since rotation is also routinely applied 
to the ultracold gases [H, , there would also be ways 
to test the scenario considered in this paper with the 
formation of vortex lattices. One could then imagine a 
mixture of clouds with bosonic particles as dark matter 
and then fermions to represent luminous matter. This 
would allow the experimental exploration of the influence 
of BEC dark matter on structures in the luminous com- 
ponent. The ability to tune self- interactions and cross 
species interaction could also allow us to test effects of 
dark matter-normal matter interactions. 

Although these experiments are extremely hard, the 
advances of recent years leave us hope that one could 
reach experimental capabilities that can explore some 
of these systems in the not too distant future. In fact, 
within traditional condensed matter the ideas of simula- 
tions for instance cosmological phenomena using exper- 
iments have been around for some [581 ] - and this type 
of thinking naturally extends to ultracold atomic gases 
which are often seen as a simulator for condensed matter 
systems. 

Alternatively, the dynamics can be simulation on a 
computer. Comparison of theory and experiment on ul- 
tracold atomic gases has shown that the Bose-Einstein 
condensates can be well described by the non-linear 
Gross-Pitaevskii equation [l8[ which can be solved nu- 
merically with great precision (see [59] for details on 
state-of-the-art numerical techniques and references to 
relevant work). The inclusion of attractive 1/r gravita- 
tional type potentials have been considered recently 
in both a variational approach and in numerical simula- 
tions. Rotation should be possible to include in similar 
fashion. We speculate that to separate dark and lumi- 
nous matter one can consider a mixtures of two bosonic 
atoms that are rotated differently so as to allow for vor- 
tices in one but not the other component or simply a 
bose-fermi mixture as also mentioned above. Issues of a 
potential thermal non-condensed part of the dark matter 
BEC could also be addressed numerically. 

A recent example of the success of the Gross-Pitaevskii 



approach is worth mentioning. The atom 52 Cr has a 
large magnetic dipole moment that can be aligned by 
externally applied fields. This means that an attractive 
dipole-dipole force can be created that was predicted to 
produce collapse of a 52 Cr condensate. This has been ex- 
perimentally observed and found to be in agreement with 
predictions from the Gross-Pitaevskii equation including 
the long-range force [6lj . This implies that long-range 
forces can be accomodated by the Gross-Pitaevskii the- 
ory as well. 



VIII. CONCLUSIONS 

We have consider some models of dark matter where 
the main component is a bosonic particle with a very 
small mass of order m ~ 10 -22 eV. Naively, this gives a 
de Broglie wavelength that is of galactic proportion such 
that one can imagine the entire galactic dark matter halo 
being in the condensate ground-state. We considered 
previous proposal with such particles with and without 
repulsive self-interactions and found general consistency 
between these for certain ranges of the mass and inter- 
action strength. Moreover, we considered the suggestion 
that superfluid BEC dark matter in rotation would likely 
also lead to vortices as seen in atomic BEC experiments. 

In case of a repulsive self-interaction wc argued that 
the vortex size should be determined locally by the co- 
herence length of this repulsive interaction. This means 
that we have two scales in the problem: A galactic one, 
given by the de Broglie wavelength from the tiny mass, 
and a sub galactic one that is determined by the mass 
and the two-body interaction strength (characterized by 
the scattering length a). We explored the consequences 
of self-interactions on the virialization of gravitationally 
bound structures and found almost no effect for reason- 
able values of m and a. 

Under the assumption of dark matter being an ultra- 
light BEC, the rotation of spiral galaxies would cause 
vortex lattices to form. We briefly addressed possible ef- 
fects on the growth of perturbations in the linear regime, 
although this is probably not relevant since the rotation 
rate in the early Universe is very small. We then consid- 
ered possible effects of sub galactic vortices in the dark 
matter on the rotation velocity curves of virialized galax- 
ies with standard dark matter halo profiles. Here we 
found that one can actually get substructure in the rota- 
tion curves that resemble some observations, but that 
this requires large vortex core size and small vortex- 
vortex distances. The mass and interaction strength 
needed to realize this where found to be fine-tuned, but 
could possibly be accommodated in more general setups. 

The present investigations and simple numerical exper- 
iments point to an interesting effect from bosonic dark 
matter. However, to fully explore the influence of vortex 
lattice formation and the feed-back on structure forma- 
tion in luminous matter one would need to consider an 
ultralight BEC dark matter component in large TV-body 
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simulations. 



Acknowledgments 

Discussions with H. 0. U. Fynbo and S. Hannestad are 
highly appreciated. Thanks also to J. P. U. Fynbo for 



[1] A. Riess et ai, Astron. J. 116, 1009 (1998). 
[2] S. Perlmutter et al, Astrophys. J. 517, 565 (1999). 
[3] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 
[4] V. V. Zhytnikov and J. M. Nester, Phys. Rev. Lett. 73, 
2950 (1994). 

[5] L. Bergstrom, Rep. Prog. Phys. 63, 793 (2000). 
[6] J. M. Overduin and P. S. Wesson, Phys. Rep. 402, 267 
(2004). 

[7] See f.x. N. J. C. Spooner, Jour. Phys. Soc. Japan 77, 

111016 (2007) and references therein. 
[8] G. Dvali, G. Gabadadze, and M. Shifman, Mod. Phys. 

Lett. A 16, 513 (2001). 

[9] R. Kallosh, [arXiv:hep-th/0205315 vl. 
[10] C. T. Hill, D. N. Schramm, and J. N. Fry, Comm. Nucl. 

Part. Phys. 19, 25 (1989). 
[11] S. J. Sin, Phys. Rev. D 50, 3650 (1994). 
[12] W. Hu, R. Barkana, and A. Gruzinov, Phys. Rev. Lett. 

85, 1158 (2000). 
[13] J. Goodman, New Astronomy 5, 103 (2000). 
[14] P. J. E. Peebles, Astrophys. J. 534, L127 (2000). 
[15] F. E. Schunck and E. W. Mielke, Class. Quantum Grav. 

20, R301 (2003). 

[16] F. E. Schunck, [arXiv:astro-ph7980 2258vl. 
[17] C. G. Bohmer and T. Harko, JCAP 06 025 (2007). 
[18] L. P. Pitaevskii and S. Stringari: Bose-Einstein Conden- 
sation, Oxford University Press Inc., New York. 
[19] M. P. Silverman and R. L. Mallett, Gen. Rel. Grav. 34, 

633 (2002). 

[20] J. Mielczarek, T. Stachowiak, and M. Szydlowski, Int. J. 
Mod. Phys.D 19, 1843 (2010). 

[21] G. Chapline and P. O. Mazur, [arXiv:gr-qc/0407033fr 2. 

[22] G. Chapline, E. Hohlfeld, R. B. Laughlin, and D. I. San- 
tiago, Int. J. Mod. Phys. A18, 3587 (2003). 

[23] E. J. Yarmchuk et ai, Phys. Rev. Lett. 43, 214 (1979). 

[24] K. W. Madison et al, Phys. Rev. Lett. 84, 806 (2000). 

[25] J. R. Abo-Shaeer et ai, Science 292, 476 (2001). 

[26] P. Engels et al, Jour. Low Temp. Phys. 134, 683 (2004). 

[27] D. O'Dell, S. Giovanazzi, G. Kurizki, and V. M. Akulin, 
Phys. Rev. Lett. 84, 5687 (2000). 

[28] C. Beauvais and G. Bothun, Astrophys. J. Suppl. 125, 
99 (1999). 

[29] I. Dicaire et al, Mon. Not. Roy. Astron. Soc. 385, 553 
(2008). 

[30] T. Matos and L. Arturo Urefia-Lopez, Phys. Rev. D 63, 

063506 (2001). 
[31] G. Baym et al, Phys. Rev. Lett. 83, 1703 (1999). 



reading the draft and suggesting valuable improvements. 



[32] A. Arbey, J. Lesgourgues, and P. Salati, Phys. Rev. D 

68, 023511 (2003). 
[33] K. Lanczos, Zeit. f. Physik 21, 73 (1924). 
[34] G. Gamov, Nature (London) 158, 549 (1946). 
[35] K. Godel, Rev. Mod. Phys. 21, 447 (1949). 
[36] S. Hawking, Mon. Not. Roy. Astron. Soc. 142, 129 

(1969). 

[37] J. Silk, Mon. Not. Roy. Astron. Soc. 147, 12 (1970). 
[38] J. D. Barrow, R. Juszkiewicz, and D. H. Sonoda, Mon. 

Not. Roy. Astron. Soc. 213, 917 (1985). 
[39] P. J. E. Peebles, Astrophys. J 155, 393 (1969). 
[40] A. G. Doroshkevich, Astrohzika 6, 581 (1970). 
[41] S. D. M. White, Astrophys. J. 286, 38 (1984). 
[42] T. Padmanabhan: Theoretical Astrophysics, Vol III, 

Cambridge University Press 2002. 
[43] M.S. Longair: Galaxy Formation, Springer- Verlag Berlin 

Heidelberg, 1998. 
[44] V. K. Tkachenko, Sov. Phys. (JETP) 22, 1282 (1966). 
[45] N. R. Cooper, E. H. Rezayi, and S. H. Simon, Phys. Rev. 

Lett. 95, 200402 (2005). 
[46] V. Schweikhard, I. Coddington, P. Engels, S. Tung, and 

E. A. Cornell, Phys. Rev. Lett. 93, 210403 (2004). 
[47] R. P. Yu and M. J. Morgan, Class. Quantum Grav. 19, 

L157 (2002). 
[48] P. Salucci, Proc. IAU Symp. 244, 53 (2007). 
[49] V. C. Rubin, D. Burstein, W. K. Ford, Jr., and N. Thon- 

nard, Astrophys. J. 289, 81 (1985). 
[50] S. M. Kent, Astron. J. 91, 1301 (1986). 
[51] A. H. Broeils, Astron. Astrophys. 256, 19 (1992). 
[52] K. G. Begeman, Astron. Astrophys. 223, 47 (1989). 
[53] J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ket- 

terle, Science 292, 476 (2001). 
[54] S.M. Carroll, Phys. Rev. Lett. 81, 3067 (1998). 
[55] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 

[56] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 

Phys. 80, 1215 (2008). 
[57] A. L. Fetter, Rev. Mod. Phys. 81, 647 (2009). 
[58] W. H. Zurek, Phys. Rep. 276, 172 (1996). 
[59] P. Muruganandam and S. K. Adhikari, Comp. Phys. 

Comm. 180, 1888 (2009). 
[60] H. Cartarius, T. Fabcic, J. Main, and G. Wunner, Phys. 

Rev. A 78, 013615 (2008). 
[61] T. Lahaye et al, Phys. Rev. Lett. 101, 080401 (2008). 



